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ABSTRACT 


This research investigates thermal buckling characteristics of rectangular panels subjected to different types of 
humped temperature profile heating. Minimum potential energy and finite-element methods arc used to calculate 
the panel buckling temperatures. The two methods give fairly close thermal buckling solutions. “Buckling 
temperature magnification factor of the first kind, if’ is established for the fixed panel edges to scale up the 
buckling solution of uniform temperature loading case to give the buckling solution of the humped temperature 
profile loading cases. Also, “buckling temperature magnification factor of the second kind, if’ is established for the 
free panel edges to scale up the buckling solution of humped temperature profile loading cases with unheated 
boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up. 

NOMENCLATURE 

Acronyms 


BC 

JLOC 

SPAR 

TPS 


A 


kl 


AR 

a 

ij 

a mnkl 

D 

D mn 

b 

r 

^ mn 


D ;; 


D Qx> °Qy 


boundary condition 
joint location 

Structural Performance and Resizing 
thermal protection system 


Symbols 

Fourier coefficient of trial function for out-of-plane displacement vv, in. 
extensional stiffness of orthotropic plates, A ^ 1 = 


A 2l = ! sVxyEy » A 22 = 7— - £ ~ Z — • A 66 = 2t s G lb/in 

1 — V V 1— VV 15 

xy yx 


t s E x A = t J y 1 £x 

1 - V *>'V 12 1 - V *.yV 


xy yx 

= a/b, plate aspect ratio 
length of plate, in. 

coefficient of characteristic equation 

Fourier coefficient of trial function for y v _ , in/in. 

width of plate, in. 

Fourier coefficient of trial function for y , in/in. 

1 yz 

EJ, 


vEJ„ 


vEJ., 


bending stiffness of plate, D ^ 1 = AA — , D 12 = — — — A_£_ ; = — A — LA 


1 — v v 

xy yx 


1 — v v 

xy yx 


1 — v v 

xy yx 


E I 


D 22 = TT7 J jr- D 6<, = 2G„/,.m-lb 

y xy yx 


transverse shear stiffness in x z-, yz-planes, lb/in. 



E 

E x , Ey 

G xy 

Is 


Young’s modulus, lb/in“ 

Young’s modulus, lb/in“ 
shear modulus, lb/in 

moment of inertia, per unit width, with respect to plate centroidal axis, / 


1 ,3 . 4/ . 
— t , in /in. 
12 5 


i 

j 

k 

I 


m 

T T T 
N N N 

x ’ y y ’ y xy 


n 

T 


T 

1 c 

T 

m 

T 

1 o 

T 

s 


v 

w 


X, y,z 


a, a x , a y 


a 


■xy 


1 xz’ Yyz 

?7 


A 

A 


c 

o 


v 


index, 1, 2, 3,.... 
index, 1, 2, 3,.... 
index, 1, 2, 3,.... 
index, 1, 2, 3,.... 

index, 1, 2, 3,...., number of buckle half waves in x-direction 
thermal forces, lb/in. 

index.l, 2, 3,...., number of buckle half waves in y-dircction 
temperature (measured from room temperature), °F 
temperature for constant temperature profile heating, °F 
material temperature 

peak temperature of dome-shaped (or roof-shaped) temperature profile heating, °F 

boundary heat sink temperature, °F 

thickness of plate, in. 

x-displacement, in. 

y-displaccment, in. 

out-of-plane displacement, in. 

rectangular Cartesian coordinates 

coefficient of thermal expansion, in/in-°F 

coefficient of thermal shear distortion, in/in-°F 

transverse shear strain in xz- and yz-planes, in/in. 

numerical factor in buckling equation written for specified panel edge condition 
= (T 0 ) ct ./ (T c ) cr ’ “buckling temperature magnification factor of the first kind” = buckling 
temperature {T Q ) cr of dome temperature loading case (fixed or free support condition) 

divided by the buckling temperature (T ) of uniform temperature loading case (fixed 

c cr 

boundaries) 

thermal buckling eigenvalue associated with constant temperature profile heating 
thermal buckling eigenvalue associated with dome-shaped (or roof-shaped) profile heating 
Poisson ratio 
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£ = [(T g ) cr ] T q /[(T 0 ) c) ] t _ Q “buckling temperature magnification factor of the second 

kind” = buckling temperature (T 0 ) of dome temperature loading case with any heat sink 
temperature T s * 0 (free support condition) divided by the buckling temperature (T g ) of 
dome temperature loading case with zero heat sink temperature T - 0 (free support 
condition) 

o 

p density, lb/in 

'y 

Oj tensile strength, lb/in 

<y 

O x stress in x-direction, lb/in 

'y 

O y yield stress, lb/in 

'y 

Oy stress in y-dircction, lb/in 

'y 

X shear stress, lb/in 

( ) cr critical value at buckling 


INTRODUCTION 

Hypersonic aircraft are subjected to severe aerodynamic heating during flights. To maintain the structural 
integrity under high temperature environment, the vehicle structural design concepts of hypersonic flight vehicles 
are different from those of low Mach number aircraft. Depending on the operating temperature range, vehicle 
structures may be called ‘“hot” structures or “warm” structures. The “hot” structures fabricated with high 
temperature alloys can operate at elevated temperatures in the range of 1000 °F to 1500 °F. If fabricated with the 
carbon/carbon composite material, the operating temperature of the “hot” structures could go as high as 3000 °F. 
Typical candidate “hot” structural components for hypersonic flight vehicles are hat-stiffened panels fabricated 
with either monolithic titanium alloys or metal-matrix composite materials; honeycomb sandwich panels 
fabricated with super-alloy; and carbon/carbon composite elevon (or body flap). 

The “warm” structures are fabricated with lightweight materials, such as aluminum, and can operate only up to 
moderate temperature limit of 350 °F. The space shuttle orbiter is a good example of the “warm” structure. The 
entire vehicle is protected with the thermal protection system (TPS) to shield the aluminum substructure from 
overheating beyond the “warm” temperature limit. 

Hot structural panels are usually fastened to the cooler substructures that function as heat sinks because of less 
heating. Thus, even under the uniform surface heat flux, the temperature distribution over the hot structural panels 
will not be uniform but looks like camel-humped shape (refs. 1-3). This camel-humped shaped temperature 
distribution is always observed in supported hot structural panels. The buckling behavior of the panel with heat 
sinks is therefore quite different from the case without heat sinks. Earlier Ko extensively studied thermal buckling 
problems of hot structural panels under uniform temperature profile heating (without heat sinks effects) 
(refs. 4-10). 

This report studies the thermal buckling behavior of rectangular panels under different types of 
humped-shaped temperature profile heating to simulate the existence of the supporting boundary heat sink and also 
studies the effect of the heat sink temperature on the panel buckling temperature. The results are compared to the 
results of uniform temperature profile heating cases without the heat sink effect. 
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HUMPED TEMPERATURE PROFILES 


Under uniform surface heating, because of the existence of the supporting cooler boundary heat sinks, 
experimental data have shown that the temperature distribution over a hot structural panel is usually 
flattened-dome shaped (fig. 1, ref. 1). For performing thermostructural analysis, the actual temperature profile may 
be idealized with simple mathematical surfaces such as dome-shaped or roof-shaped temperature profiles 
described in the following sections. 

Dome-Shaped Temperature Profile 

The dome-shaped temperature profile for a rectangular panel (length a, width b, thickness t s ) is described by 
the following mathematical function: 


T(x,y) = T +{T - T )sin — sin^ 
^ os a b 


( 1 ) 


where T a is the peak temperature at the panel center and T s is the boundary heat sink temperature (fig. 2). 
Equation (1) gives a “sine-sine” surface elevated upward by an amount T s (fig. 2). Keep in mind that the 
temperatures in this report are the temperature differentials above the room temperature. 

Equation (1) is the first-order approximation of the actual temperature profile (fig. 1). For more accurate 
representation of the actual temperature profile, double Fourier series representation is necessary; however, this 
representation is not the scope of the present report. 

When the boundary temperature T s goes to zero (fig. 3), equation (1) becomes: 

T(x,y) = T 0 sin^sin^ (2) 

a b 

In order to obtain closed-form buckling solutions by means of the minimum potential energy theory, the functional 
form given by equation (2) is used for dome-shaped temperature loading. 

For a special case of uniform temperature profile heating (fig. 4), equation (1) degenerates into: 

T(x,y) = T (> = T s = T c (3) 

Thermal buckling problems of rectangular sandwich panels under the constant temperature profile heating 
[eq. (3), fig. 4] was analyzed earlier by Ko (refs. 4—7) using the minimum potential energy theory and 
finite-element method. 


Roof-Shaped Temperature Profile 

The actual dome-shaped temperature profile (fig. 1) may also be approximated by a simple temperature 
distribution surface called roof-shaped temperature profile (peak temperature T a at the profile central flat zone and 
boundary heat sink temperature T s ) (figs. 5 and 6). When T Q = T , the roof-shaped temperature profile also 
degenerates into uniform temperature profile (fig. 4). 
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PROBLEMS INVESTIGATED 


This report investigated thermal buckling behavior of rectangular panels (length a, width b, thickness t s ) 
heated under the dome-shaped (figs. 2 and 3) and roof-shaped (figs. 5 and 6) temperature profiles. The main 
aspects of this research are: 

1 . To study the effect of the dome-shaped (or roof-shaped) temperature profile on the buckling behavior of the 
rectangular panels with fixed edges. Then, the study will establish the “buckling temperature magnification 
factor of the first kind” to scale up the buckling solution of uniform temperature profile heating case to give 
the buckling solution of the dome -/roof-shaped temperature profile heating case. 

2. To study the effect of the heat sink temperatures on the buckling temperatures of rectangular panels with 
free edges under dome-shaped (or roof-shaped) temperature profile heating. Then, the study will establish 
the “buckling temperature magnification factor of the second kind” to scale up the buckling solution of 
dome (or roof) temperature loading case with unheated boundary heat sinks to give the buckling solutions 
when the boundary heat sinks are heated up. 


FINITE-ELEMENT BUCKLING ANALYSIS 


This research will use the Structural Performance and Resizing (SPAR) finite-element computer program (ref. 
1 1) to conduct the 1 i near elastic thermal buckling analysis of the rectangular panels. 

Panel Geometry 

All the rectangular panels analyzed have thickness t s - 0.09 in. with the different dimensions in length a and 
width b listed in Table 1. 


Table 1. 

Dimensions of rectangular panels; t s 

= 0.09 in. 

a, in. b, in. 

alb 

6 

6 

1.0 

9 

6 

1.5 

12 

6 

2.0 


Figure 7 shows the quarter model generated for the rectangular panel. The quarter model has 625 joint 
locations (JLOCs) and 576 E43 four-node plate elements. 

Material Properties 

The rectangular panels are made of high temperature alloy Flaynes® 230® (Flaynes International, Koyomo, 
Indiana) with the room temperature properties listed in Table 2. 
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Table 2. Material properties of Haynes 230 alloy at room temperature. 


E 

30.6 x 10 6 lb/in 2 

V 

0.31 

P 

0.324 lb/in 3 

a 

7.0 x 10~ 6 in/in- °F 


Haynes 230 high temperature alloy was used in the fabrication of Hyper-X (X-43) hypersonic research vehicle 
for Mach 7-10 mission (ref. 3). Appendix A gives the temperature-dependent material properties of Haynes 230 
alloy. 


Boundary Conditions 

The following different panel boundary conditions were considered in the finite-element analysis: 

1. 4S fixed — four edges simply supported, with no in-plane displacements [fig. 8(a)]. 

a. u(x,y) ^ 0, v(x,y) ^ 0 everywhere except boundaries. 

b. u(x,y) = 0, v(x, y) = 0 everywhere in the panel. 

2. 4S free — four edges simply supported and slides freely along the lubricated guides, which can have free 
in-plane motions [fig. 8(b)]. 

3. 4C fixed — four edges clamped, with no in-plane displacements [fig. 9(a)]. 

a. u(x,y) ^ 0, v(x,y) ^ 0 everywhere except boundaries. 

b. u(x,y ) = 0, v(x, y) = 0 everywhere in the panel. 

4. 4C free — four edges clamped and slides freely along the lubricated clamping guides, which can have free 
in-plane motions [fig. 9(b)]. 


Thermal Loads 

The thermal load inputs used for the dome (or roof) temperature loading and constant temperature loading are 
as follows. 

1 . Dome-Shaped (or Roof-Shaped) Temperature Profile Heating 

For the dome-shaped [eqs. (1) and (2), figs. 2 and 3] and roof-shaped (figs. 5 and 6) temperature load inputs to 
the SPAR program for the eigenvalue calculations, the unit peak temperature T () = 1 °F was used. The heat sink 
temperature T s was allowed to vary over the range T S /T ( = 0~1 (i.e., T g /T Q = 0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 
0.8, 0.9, 1.0). Then the eigenvalue (scaling factor) X () calculated from SPAR program will give the buckling 
temperature (T Q ) , namely: 


(T ) = A x 

v o'er o 


1 = A 


(4) 
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If the input peak temperature T a is doubled (i.e., T () = 2 °F), then the calculated eigenvalue will be reduced to half 
of A given in equation (4) because the product (T x A ) remains constant. 

The reason for using the whole range of T /T = 0~1 is to study the effect of the heat sink temperature T s on 
the panel buckling temperature (T ) when the panel edges are free to move. Keep in mind that the panel 
buckling temperature (T 0 ) will increase with the increasing heat sink temperature T s because of relaxation of 
thermal expansion constraint exerted on the panel by the boundary heat sink. 

2. Uniform Temperature Profile Heating 

As mentioned earlier, the uniform-temperature loading case (fig. 4) is a degenerative case of dome-shaped (or 
roof-shaped) temperature profile heating. Temperature load of T c = 1 °F is chosen as input to all nodes of the 
finite-element model so that the eigenvalue A f . calculated from SPAR program will give the buckling temperature 

( r c) 0 -’ namely: 


Wc) cr = K*' = K ( 5 ) 

This buckling solution is for the panels with fixed boundary supports. When the panel boundaries can have free 
in-plane motions (free supports), the constant temperature loading case obviously can never induce thermal 
buckling (i.e., buckling temperature goes to infinity). 

3. Material Property Iterations 

In the calculations of buckling temperatures, material property iterations are required to obtain accurate 
buckling temperatures. Namely, the material properties at certain assumed material temperature T m must be 
updated until T m matches the calculated buckling temperature (T c ) . Reference 12 discusses in detail this 
iteration process. Usually it takes only two to three material iterations to yield accurate buckling temperatures. 

BUCKLING TEMPERATURE MAGNIFICATION FACTORS 


As discussed later in this report, for the fixed supports (4S or 4C) [figs. 8(a) or 9(a)], the buckling temperature 
(T g ) cr for the dome (or roof) temperature loading case is much higher than the buckling temperature (T c ) of 
the uniform temperature loading case. Therefore, a “buckling temperature magnification factor of the first kind, rf’ 
defined as: 


V = 


( T ) 

v o'er 

(T r ) 
c' cr 


(6) 


will be used to indicate how many times the buckling temperature (T ) of the dome (or roof) temperature 

loading case (fixed supports) is magnified from the buckling temperature (T ) of the uniform temperature 

c cr 

loading case with fixed supports (fundamental case). 

In reality, the hot panel attempts to expand under heating; but its expansion is resisted by the cooler boundary 
substructures (heat sinks) that expand less. This boundary constraint is the cause of thermal buckling of the panels. 
Such substructure constraints will gradually relax as the substructures are heated up, resulting in higher panel 
buckling temperature. In order to discover how the heat sink temperature T s affects the panel buckling temperature 
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(T 0 ) for the dome (or roof) temperature loading case with free supports, another “buckling temperature 
magnification factor of the second kind, defined as: 


(T ) 
v o’er 


0 


(T ) 

V o’er 


Js = 0 


(7) 


will be used to indicate how many times the buckling temperature {T 0 ) cr for the dome (or roof) temperature 
loading case with free support (4S or 4C) is magnified when the heat sink temperature T s increases from T = 0 
(no heat sink thermal expansion) to certain nonzero value T * 0 (with heat sink thermal expansion). 


MINIMUM POTENTIAL ENERGY BUCKLING ANALYSIS 


When the {«, v} displacements are zero everywhere in the plate including the boundaries (4S fixed-b and 4C 
fixed-b cases), the dome-shaped temperature profile heating [eq. (2), fig. 3] will generate identical dome-shaped 
biaxial stress distributions with no shear stresses. Because of this ideal condition, it is possible to obtain 
closed-form buckling solutions for this particular case using the minimum potential energy theory. 

Energy Equations 

Strain energy associated with bending V\ may be written as (refs. 5 and 6): 


Vl = 


a i 

a 


b 1 D r 


d / dw _ 
dx\ dx ^ xz > 


-|2 


+ D 


12 


d / dw _ 




Vdx V dx xz ) ] Ldy \ dy yz ) J 


d tdw _ 




D 21[ d (dw 


+ ~ L^"MJ + 2 


-|2 D 


66 


A(^-y \+±(*Z-y 

.dy\dx xz ) dx\ dy yz 


-|2 D 


Qly 2 

2 xz 


D 


+ —Y X7 +^ 1 Y V7 \dxdy 


2>y 2 1 

2 ryz j 


(8) 


and the strain energy associated with thermal loading V 2 is given by: 

a b r t / j,.,\ 2 




(9) 


The dome-shaped temperature loading functions may be expressed as follows (refs. 5 and 6): 


N x (x,y) = (A n a x + A l2 a y )T o sin — sin 


jcx ■ icy 

cm 

b 


( 10 ) 


N v (x,y) = (A 2 ia x + A 2 2Cc y )T 0 sm—sm^ 


b 


( 11 ) 


N xy (x,y) = 0 


( 12 ) 


The preceding thermal loading functions (10)— (12) hold only for the 4S fixed-b and 4C fixed-b cases 
[i.e., u(x,y) = 0, v{x,y) = 0 everywhere]. 
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For other fixed panel cases under dome temperature heating, the in-planc displacements of panel interior points 
are nonzero [i.e., u(x, y) * 0, v(x,y) * 0], and the thermal loads { N x , N y , N xy } are complicated functions of {x,y}\ 
and thus it will be too cumbersome to apply the minimum potential energy method. Therefore, SPAR was used to 
obtain quick buckling solutions. 


Deformation Functions 

The deformation functions for the 4S fixed-b case and 4C fixed-b case may be expressed in the following 
forms (refs. 5 and 6). 

Case 1. Four edges simply supported (4S fixed-b case) 


w{x,y) 


CO 00 


V' v'' i ■ mnx ■ nny 

\ \ 4 . sin sin — *- 

L Z mn a b 

m = 1 n = 1 


GO 00 


/ \ v'' V'' /> mnx ■ nny 

y X z^ = 2 Z B 'nn COS — Sm 


m — \ n — 1 


00 00 


= 2 2 c 


■ mux mty 

sin cos — - 

mn a b 


m = In = 1 


Case 2. Four edges clamped (4C fixed-b case) 


w{x,y) = sin ^ sin & J J A m« sin 


m xx ■ nny 

' sin — 

a b 


1 n = 1 


00 00 

/ , nx ■ ny s-' v* i, ■ mnx ■ nny 

rjx.y'i = S "'T 

m = 1 n = 1 

00 00 


, ■ nx ■ ny 

+ sin — sin — 
a Z? 


\-\ v' ,, mnx ■ nny 

\ \ mo,,,, c os sin — i 

Z Z a b 

m = lw = 1 


00 CO 


/ \ tdc Try v' v' ^ 

y yz (*,)9 = sm ~ cos_ ^ 2j 2j C 


m = 1 n = 1 

00 00 


, ■ nx ■ ny 

+ sin — sin — 

a b 


■ mnx ■ nny 

sin sin — • L 

a b 


mnx nny 

\ nr\ c 4 


V' fx nuijy r 

\ \ nC „ sin cos- 

Z Z mn a b 


m = In = 1 


(13) 


(14) 


(15) 


(16) 


(17) 


(18) 


Buckling Equations 

Application of the minimum potential energy theory yields the following homogeneous simultaneous 
characteristic equations (refs. 5 and 6) for each set of integral values {m, n } (or mode shape) written for orthotropic 
rectangular panels under dome-shaped temperature profile heating [eq. (2)]: 
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I? 


M 


mnkl 


+ P 


mnkl 


A k , = 0 


(19) 


where the stiffness/geometry parameter M mnk[ is given by: 


M 


mnkl 


12 . 23 31 21 33 . 13 , 21 32 22 31 ' 

11 , a mnkA a mnkl a mnkl a mnkl a mnkV a mnkl' a mnkl a mnkl a mnkl a mnkV 


a mnkl 


22 33 23 32 

a mnkl a mnkl a mnkl a mnkl 


(20) 


and where a J , , are the coefficients of characteristic equations. 


Appendix B gives functional expressions of . , , which is taken from Appendix C of reference 6 (or 
Appendix A of ref. 7) with the loading terms containing {k x , k Y ] removed. 

For the simply-supported edges (4S fixed-b case), the value of 'Q in equation (20) is £ = 32; and the thermal 
forcing term in equation (19) has the following form (special integrals given in Appendix C were used in the 
calculations). 


mnkl 


'j j Am 2 + k 2 )- 1 — An+l 2 )- 1 

(A ll a x + A n a y ) + (^2l a x + A 22 a y ) 3 

a b 


( 21 ) 


mnkl 


l(m + k) 2 - l][(m-k) 2 -!][(« + I) 2 - !][(«- if- 1], 


m ± k - even , n ± l = even 

For the clamped edges (4C fixed-b case), the value of l, in equation (20) is £= 24, and in equation (19) 
takes on the following form (special integrals given in Appendix C were used in the calculations). 


_ mnkl 
mnkl ~ 2 

a 


(^ll a x + ^12 a y) 


3 {m 2 + k 2 )+ 1 


3 [(m +k“)-5] 


■ V \(m + k) 2 -l][(m-k) 2 -l] [(m + k) 2 - 9][(m - k) 2 - 9] . 


( 22 ) 


-1 


{[(n + l) 2 - 1 ][(n - if - 1 ] [(» + if -9 ][(« - if - 9] . 


3(« 2 + / 2 )+ 1 


3[(n 2 + / 2 )-5] 


, mnkl , . , , \ I 

+ 2 ( A 21 a x + A 22 a yf ? 2 2 2 

b {[(n + l) -1] [(n + /) -9][(n-/)“-9], 


-1 


+ 


\(m + kf - 1 ][(m - kf - 1 ] \{m + k) -9][(m-k) -9], 


m ± k - even , n ± 1 = even 
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For the unknown deflection coefficients in equation (19) to have nontrivial solution, the determinant of the 
coefficients of unknown A k j of simultaneous equations written out from equation (19) must be set to zero for 
eigenvalue extractions. The largest eigenvalue (1 IT 0 ) will then give the lowest critical buckling temperature (T 0 ) cr . 

Appendix D shows the determinants of the coefficients of simultaneous equation (called buckling equations) 
written out from equation (19) up to order 12 for the cases m + n = even (symmetrical buckling) and m ± n = odd 
(antisymmetrical buckling) for the 4S and 4C cases. The symmetrical buckling case (m ± n = even) gives the 
lowest buckling temperature (T 0 ) cr associated with the lowest buckling mode. 

Buckling Temperatures for Isotropic Cases 

The purpose of the present buckling analysis is to study the effect of dome-shaped temperature profile heating 
on the buckling of the rectangular panels and not the effect of material anisotropy. Therefore, only the isotropic 
panels are considered. 

1. Simply-Supported Edges 

For the simply-supported case (4S fixed-b case) ignoring the transverse shear effect, M mrM [eq. (20)] and P mn ki 
[eq. (21)] take on the following forms: 


M 


T7? 4 

Et n 


mnkl 9 

384(1 -v ) 


2 \ 


2 

m n 

— + 72 

a b 


(23) 


mnkl 


Et $ a 

1 - V 


(m + k~) - 1 (n~ + r)— 1 


(24) 


mnkl 

X{m + k) 2 - l][(m-k) 2 -!][(« + if- !][(«- if- 1], 


m±k = even, n + l = even 

The first-order solution (m = n = k = l = 1) of the buckling temperature (T 0 ) cr for the dome-shaped heating may 
be calculated through combining equations (19), (23), and (24) by setting m = n = k= l = 1. 


v °’cr 128 a(l 


(25) 


For the constant temperature heating case, the first-order (in = 1, n = 1) buckling temperature (T c ) cr for the 
isotropic panel can be calculated from equations (28)-(30) of reference 7. With the transverse shear effect 
neglected, the first-order buckling temperature (r c ) c ,.for this case has the following form: 


11 



(26) 


2 2 
It t 

( T ) = 

c cr 12 a( 1 + v) 


1 + 1 

2 . 2 
a b 


\ 


Dividing equation (25) by equation (26) yields the buckling temperature magnification factor of the first kind 
?7 for the simply- supported first-order case (in -n-k = l= 1) in a neat form: 


V = 


( T ) 

v o'er 

( r c ) 

c cr 


9 It 
32 


= 2.7758 


(27) 


that is independent of panel aspect ratio. 

For the 4S case, the values of ?7 based on the second and third-order buckling solutions (which depend on the 
aspect ratio alb ) were also obtained. The RESULTS section presents these values. 

2. Clamped Edges 

For the clamped case (4C fixed-b case), the first-order solution (m = n - k = 1 = 1) for the “dome-shaped” 
temperature case is: 


M 


4 3 

it Et , 


till 


24(1- v ) 


! + 2_L + ± 

4 3 2,2 ,4 

a abb 


\ 


(28) 


512 


Eat,, 


till 


675(1 -v) 


\ 


1 + 1 

a 2 b 2 


(29) 


The first-order solution (m = n = k = 1 = 1) of the buckling temperature (T 0 ) cr for the clamped panel under the 
dome-shaped heating may be calculated by combining equations (19), (28), and (29) to yield: 


(T ) = 

v O’r.r 


4 2 

15nt s 


^3 2 3 

— + — — + — 

4 2 2 4 

\a a b b 


64 ~a( 1 + v) 


T + l) 

\« 2 b 2 j 


(30) 


The first-order solution (m = n - k = l - 1) for the constant temperature case calculated from equations 
(28)-(30) of reference 7 gives: 


2 2 
It t „ 

(T ) = — 

c cr 9a(l + v) 


^3 2 3 

— + + — 

4 ii a 
\a a~b b 


T + J-i 

\a 2 b 1 ) 


(31) 
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The ratio of equations (30) and (31) gives the buckling temperature magnification factor of the first kind ?7 for 
the 4C case first-order solution (m = n = k= 1 = 1) in a neat closed form as: 

(T 0 ) 3 x i ^ 2 ^ 2 

t] = —^1 = - x 13 - = 1.6265 (32) 

( T d r 64 2 

which, similar to the 4S case [eq. (27)], is also independent of panel aspect ratio. As discussed in a following 
section, for the 4C case, the first-order buckling solution gives a quite accurate value of //. Therefore, higher order 
values of ?7 were not calculated. 

RESULTS 

The following sections present the results of the thermal buckling analyses of the rectangular plates under 
different heating profiles. 

Buckling Temperature Magnification Factors 

Figures 10 and 1 1 show the buckled shapes of the isotropic square panel ( alb =1) under different temperature 
profile loading and under different boundary conditions. Notice that the panel buckling shapes are insensitive to 
the temperature loading functions for each set of edge conditions. 

Table 3 lists the buckling temperature magnification factor of the first kind // = () ) (] / {T c ) (r calculated 

from SPAR program for the isotropic rectangular panels with different aspect ratios under dome-shaped 
temperature profile heating with zero heat sink temperature \T S - 0, eq. (2)]. 

Table 3. Buckling temperature magnification factor ?7 
for rectangular plates under dome-shaped temperature 
profile heating; T s = 0; fixed boundaries. 




?7 


BC \ alb 

1 

1.5 

2 

4S 

1.9721 

(2.0034) 

1.9603 

1.9453 

4C 

1.8035 

1.7848 

1.7320 


(1.8299) 




( ) Sandwich panel (ref. 7). 

Notice that ?7 is insensitive to the change of panel aspect ratio alb and that the 4S cases have slightly higher 
values of ?7 than the 4C cases. 

In Table 3, the ?7 values for the square sandwich panel (alb = 1) are also shown in the parentheses. Those 
values were calculated from the previous SPAR program developed for the buckling of sandwich panel (ref. 7). 
Notice the proximity of the ?7 values for solid and sandwich panels. 
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Table 4 compares the values of ?7 calculated from the minimum potential energy theory and from the SPAR 
program for the special case of square panel {alb =1) under the dome-shaped temperature profile heating \T S = 0, 
eq. (2), fig. 3], and with constraint u - v = 0 everywhere. 


Table 4. Buckling temperature magnification factor ?7 for square plate ( alb = 1) 
under dome-shaped temperature profile heating; T s = 0; u = v = 0 everywhere. 


BC Method 


?7 


Ko theory 

SPAR 

Percent difference 

4S 

2.4503*** 

2.5185 

2.71 


(2.6394)** 


4.80 


(2.7758)* 


10.22 

4C 

1.6265* 

1.5993 

1.70 


* First-order solution. 

** Second-order solution. 
*** Third-order solution. 


For the 4S case, the ?7 values were hand calculated up to the third-order buckling solutions. The third-order 
solution compares quite well with the SPAR solution with only 2.71 percent difference. Thus, higher order 
solutions were not pursued. 

For the 4C case, the first-order solution is extremely good with only 1.70 percent difference compared with the 
SPAR solution. Therefore, higher order solutions for ?7 were not calculated. 

Table 5 lists the buckling temperatures {T 0 ) cr and the values of buckling temperature magnification factors 
{? 7 > s! calculated for the simply-supported (4S fixed/free) square panel ( alb =1) under shifted dome-shaped 
temperature loading [eq. (1), fig. 2] under different heat sink temperatures T s . 
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Table 5. Buckling temperatures of square plate ( alb = 1) under 
dome-shaped temperature profile heating; 4S fixed and free cases. 


T IT 

± s' 

4S Fixed 

4S Free 


(T ) °F 

v-* o'er ’ 

?7 

(T ) °F 

o’er ’ 

s 

0.0 

80 

2.00 

515 

1.00 

0.1 

73 

1.83 

568 

1.10 

0.2 

67 

1.68 

635 

1.23 

0.3 

62 

1.55 

718 

1.39 

0.4 

57 

1.43 

821 

1.59 

0.5 

54 

1.35 

960 

1.86 

0.6 

50 

1.25 

1,175 

2.28 

0.7 

47 

1.18 

1,500 

2.91 

0.8 

45 

1.13 

2,070 

4.02 

0.9 

42 

1.05 

(4,213) 

8.18 

1.0* 

40 

1.00 

00 

00 


Uniform temperature loading. 

( ) Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Note that for the 4S fixed edges, the value of )] decreases with the increasing heat sink temperature T s . For the 
4S free edge case, £ increases with the increasing heat sink temperature T s . When TJT a = 1 , the case degenerated 
into the uniform temperature loading case for which the buckling temperature goes to infinity for the 4S free edge 
case (i.e., thermal buckling does not occur). 

It must be mentioned that for the actual hot sandwich panel (ref. 1) TJT a = 0.45 and for the FIYPER-X wing 
panel (ref. 3) TJT 0 = 0.54. 

Table 6 lists the buckling temperatures (T 0 ) cr and the values of buckling temperature magnification factors 
{? 7 , s! calculated for the clamped (4C fixed/free) square panel {alb =1) under shifted dome-shaped temperature 
loading [eq. (1), fig. 2] and under different heat sink temperatures T s . 
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Table 6. Buckling temperatures of square plate ( alb = 1) under 
dome-shaped temperature profile heating; 4C fixed and free cases. 


T IT 

± s' ± o 

4C Fixed 

4C Free 


(T ) °F 

?7 

(T ) °F 

s 

0.0 

194 

1.81 

810 

1.00 

0.1 

180 

1.68 

886 

1.09 

0.2 

167 

1.56 

982 

1.20 

0.3 

156 

1.46 

1,110 

1.34 

0.4 

147 

1.37 

1,272 

1.57 

0.5 

138 

1.29 

1,480 

1.83 

0.6 

130 

1.21 

1,770 

2.19 

0.7 

124 

1.16 

2,302 

2.84 

0.8 

118 

1.10 

(3,581) 

4.42 

0.9 

112 

1.05 

(6,907) 

8.53 

1.0* 

107 

1.00 

00 

00 


Uniform temperature profile heating. 

( ) Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Note that for the 4C fixed case, similar to the 4S fixed case, value of ?7 decreases with the increasing heat sink 
temperature T a . For the 4C free case, the values of £ are very close to the corresponding £ value of 4S free case 
(Table 5). Again, when T s reaches T a (i.e., TJT a = 1), the case degenerated into the uniform temperature loading 
case for which the buckling temperature goes to infinity for the 4C free case (i.e., no thermal buckling). 

Table 7 lists the buckling temperatures (T 0 ) cr and the values of buckling temperature magnification factors 
{? 7 , q} calculated for the simply- supported (4S fixed/free) square panel {alb = 1) under shifted roof-shaped 
temperature loading (fig. 5) under different heat sink temperatures T s . 
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Table 7. Buckling temperatures of square plate (alb =1) under 
roof-shaped temperature profile heating; 4S fixed and free cases. 


T IT 

* s' ± o 

4S Fixed 

4S Free 


(T ) °F 

?7 

(T ) °F 

£ 

0.0 

51 

1.28 

622 

1.00 

0.1 

50 

1.25 

690 

1.11 

0.2 

48 

1.20 

767 

1.23 

0.3 

47 

1.18 

860 

1.38 

0.4 

46 

1.15 

985 

1.58 

0.5 

45 

1.13 

1,160 

1.86 

0.6 

44 

1.10 

1,410 

2.27 

0.7 

43 

1.08 

1,764 

2.84 

0.8 

42 

1.05 

(2,587) 

4.16 

0.9 

41 

1.03 

(5,174) 

8.32 

1.0* 

40 

1.00 

00 

00 


Uniform temperature profile heating. 

( ) Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 


Table 8 lists the buckling temperatures (T 0 ) cr and the values of buckling temperature magnification factors 
{? 7 , s} calculated for the clamped (4C fixed/free) square panel ( alb =1) under shifted roof-shaped temperature 
loading (fig. 5) under different heat sink temperatures T s . 


Table 8. Buckling temperatures of square plate ( alb =1) under 
roof-shaped temperature profile heating; 4C fixed and free cases. 


TJT a 

4C Fixed 

4C Free 

(T ) °F 

?7 

(T ) °F 

£ 

0.0 

132 

1.23 

1,125 

1.00 

0.1 

129 

1.21 

1,235 

1.10 

0.2 

126 

1.18 

1,370 

1.22 

0.3 

123 

1.15 

1,520 

1.35 

0.4 

121 

1.13 

1,720 

1.53 

0.5 

118 

1.10 

1,980 

1.76 

0.6 

116 

1.08 

(2,501) 

2.22 

0.7 

113 

1.06 

(3,335) 

2.96 

0.8 

111 

1.04 

(5,002) 

4.45 

0.9 

109 

1.02 

(10,004) 

8.89 

1.0* 

107 

1.00 

00 

00 


Uniform temperature loading. 


( ) Exceeded melting range 2,375 °F~2,500 °F of Haynes 230 alloy. 
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Figures 12 and 13, respectively, show the buckling temperatures (T 0 ) cr (data from tables 5-8) plotted as 
functions of normalized heat sink temperature TJT 0 for the fixed and free edge cases to illustrate entirely different 
buckling behavior of the fixed and free edge cases. Note from figure 12 that for the fixed edges, the buckling 
temperatures of the roof-shaped temperature profile case are much lower than the dome-shaped temperature profile 
case for each panel support condition. Flowever, for the free edges (fig. 13), the reverse is true. 

Figure 14 shows the buckling temperature magnification factors { 77 . £} (data from tables 5-8) plotted as 
functions of normalized heat sink temperature TJT a for easy visualization of the shapes of { 17 , £} curves. Note that 
?7 values for the roof case are much lower than the dome case for the same support condition. For the free edge 
cases, the £ curves for different loading and edge support conditions stay very close. 

In-Plane Deformations 

When the panel is constrained only at the four boundaries, the panel interior points can have in-plane free 
motions if the temperature loading is nonuniform. Figure 15 shows the in-plane deformed shapes of the square 
panel under dome-shaped [fig. 15(a)] and roof-shaped [fig. 15(b)] temperature loading cases. As expected, the 
panel center regions expand more than the panel boundary regions as the result of heating temperature profiles. 
Figures 16 and 17, respectively, show the distributions of in-plane displacements u(x,y) and v(x,y) for the 
dome-shaped temperature loading. Note that { 11 , v] reach a maximum at points between the panel center and the 
boundaries. For the roof-shaped temperature profile case (figs. 18 and 19), the maximum points of {u, v] migrated 
toward the panel center. 


Thermal Stresses 

Figures 20 and 2 1 , respectively, show the compressive thermal stresses { o x , o } induced in a fixed square 
panel (alb =1) under the unit dome-shaped temperature profile heating \T a = 1 °F, T s = 0; eq. (2), fig. 3], Even 
though the loading temperature profile is dome shaped, the distributions of thermal stresses { 0 , 0 } are 

x y 

airplane-hanger shaped because u(x,y) ^ 0 and v(x,y)^0 within the fixed boundaries. Figure 22 shows the 
distribution of shear stress T xy which is wavy shaped. The peak magnitudes of the shear stress occur at the 
boundaries and near the panel corners. 

For the roof-shaped temperature profile heating (T 0 = 1 °F, T s = 0) (fig. 6), the distributions of { o x , o y } are 
roughly trapezoidal hanger-shaped (figs. 23 and 24). The shear stress distribution (fig. 25) exhibits zero shear 
stress near the panel center region and a maximum shear near each panel corner. 

If the in-plane displacements {u, v] are constrained everywhere in the panel including the edges 

[i.e., u(x,y) = v(x,y) = 0], the resulting distributions of the biaxial stresses { cr , 0 } will be dome-shaped [ust like 

x y 

the input dome temperature distribution (figs. 26 and 27), and the in-plane shear stress t diminishes ( r vv = 0). For 
the case of roof-shaped temperature profile (figs. 28 and 29), the distributions of the biaxial stresses { o x , o y } 
reflect the input temperature profile. 

Table 9 lists the fundamental thermal stresses induced in the fixed square panel ( alb = 1) by the dome-shaped 
(or roof-shaped) temperature loading (T a = 1 °F) with zero heat sink temperature (T s = 0). 
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Table 9. Peak thermal stresses { o x , (7 y , T xy } induced in square plate (alb =1) under 
dome-shaped and roof-shaped temperature loads; T a = 1 °F; T s = 0; fixed boundaries. 


Constraints of {u, v} 

Ideating 

O x , lb/in 2 

O y lb/in 2 

T xy lb/in 2 

u = 0, v - 0, edges only 

Dome 

Roof 

-189.26 

-256.38 

-189.26 

-256.38 

31.27 

48.66 

u = 0, v - 0, everywhere 

Dome 

Roof 

-301.05 

-306.77 

-301.05 

-306.77 

0.00 

0.00 


Note that by setting u(x,y) = v(x,y) = 0 everywhere in the panel including the boundaries, the peak magnitudes 
of the biaxial stresses { o x , cr v } were raised respectively by 59 percent and 20 percent for the dome- and 
roof-shaped temperature profile cases. 


CONCLUDING REMARKS 

Thermal buckling characteristics of rectangular panels subjected to dome-shaped and roof-shaped temperature 
profile heating were investigated using the finite-element method and the minimum potential energy method. The 
key results include: 

1. “Buckling temperature magnification factor of the first kind, ? 7 ” was established to scale up the buckling 
solution of uniform temperature loading case to give the buckling solution of the dome-shaped (or 
roof-shaped) temperature loading cases. Also, “buckling temperature magnification factor of the second 
kind, was established to scale up the buckling solution of dome temperature loading case with unheated 
boundary heat sinks to give the buckling solutions when the boundary heat sinks are heated up. 

2. For the fixed boundary cases when the panel interior in-plane motions are not constrained (u * 0 , v ^ 0 
except boundaries), the panel buckling temperatures under dome-shaped temperature loading are 
practically twice the buckling temperatures of the uniform temperature loading cases. 

3. For simply-supported case (4S fixed, u - v = 0 everywhere) under the dome-shaped temperature profile 
heating, the third-order solution of the “buckling temperature magnification factor of the first kind, > 7 ” 
calculated from the minimum potential energy theory, agrees fairly well with the ?7 value calculated from 
the finite-element method. The solution difference is only 2.71 percent. 

4. For the clamped case (4C fixed, u = v = 0 everywhere) under dome-shaped temperature profile heating, the 
first-order solution of minimum potential energy gives quite accurate value of ?7 with only 1.70 percent 
difference from the finite-element solution. 

5. For the fixed boundary case when the panel interior points can have free in-plane motions (u * 0 , v ^ 0 
except boundaries), the distribution profiles of the compressive stresses { o x , cr y } under the dome-shaped 
(or roof-shaped) temperature loading have airplane-hanger-like shapes. 

6 . For the fixed boundary case when the panel in-plane motions are constrained everywhere (u = v - 0 
everywhere), the distribution profiles of the compressive stresses { o x , cr v } take on the dome shape (or roof 
shape) similar to the input dome-shaped (or roof-shaped) temperature profile. 

7. The “buckling temperature magnification factor of the first kind, 77 ” associated with the dome-shaped 
temperature profile heating is much greater than that for the roof-shaped temperature profile heating case 
for both 4S and 4C fixed edge conditions. 
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8. The “buckling temperature magnification factor of the second kind, associated with the dome-shaped 
and roof-shaped temperature profile heating are very close for 4S and 4C free edge conditions. 


Dryden Flight Research Center 

National Aeronautics and Space Administration 

Edwards, California, April 3, 2002 
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APPENDIX A 

TEMPERATURE-DEPENDENT MATERIAL PROPERTIES OF 

HAYNES 230 ALLOY 


o 

Temperature-dependent material properties of Haynes 230; p = 0.324 lb/in . 


T, °F 

Ex 10 6 , lb/in 2 

o T x 1 ()’ , lb/in 2 

Oy x 10 3 , lb/in 2 

a x 10 6 , in/in-°F 

* 

V 

70 

30.6 

125.4 

57.4 

7.0 

0.310 

200 

30.1 

(122.2) 

(55.0) 

7.1 

0.311 

400 

29.3 

(117.3) 

(51.3) 

7.2 

0.315 

600 

28.3 

(112.3) 

(47.7) 

7.4 

0.318 

800 

27.3 

(107.4) 

(44.0) 

7.6 

0.321 

1000 

26.4 

102.5 

40.3 

7.9 

0.324 

1200 

25.3 

97.7 

39.5 

8.1 

0.330 

1400 

24.1 

87.7 

42.5 

8.3 

0.332 

1600 

23.1 

63.1 

37.3 

8.6 

0.334 

1800 

21.9 

35.2 

21.1 

8.9 

0.340 

2000 

20.7 

19.5 

10.8 

* 

9.2 

0.343 


Estimated. 


( ) Interpolated. 
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APPENDIX B 

COEFFICIENTS OF CHARACTERISTIC EQUATIONS 


The characteristic coefficients aj r , , appearing in equation (20) are defined in the following for different 
indicial and edge conditions (refs. 6. 7). 

Case 1: 4S condition 

1. m =k,n =1 
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APPENDIX C 

TABLE OF SPECIAL INTEGRALS 


Different types of integrals are needed to carry out the integration of strain energy expression associated with 
the dome-shaped temperature profile heating. 
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APPENDIX D 
BUCKLING EQUATIONS 
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Case 1: 4S condition 
m ± n = odd (antisymmetric buckling) 
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Case 2: 4C condition 


m + n - even (symmetric buckling) 
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Case 2: 4C condition 
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FIGURES 



Figure 1. Measured temperature distribution in upper surface of titanium honeycomb-core sandwich panel, heated 
on upper side at 10 °F/sec heating rate, with four edges supported by test fixtures (heat sink) (ref. 1). 



Figure 2. Shifted dome-shaped temperature profile heating ( T * 0). 
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Figure 8. Two types of boundary conditions for simply-supported edges. 
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Figure 9. Two types of boundary conditions for clamped edges. 
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(a) Uniform temperature profile heating. 



(b) Dome-shaped temperature profile heating. 



(c) Roof-shaped temperature profile heating. 


Figure 10. Buckled shapes of square panel {a/b = 1) under different temperature profile heating; 4S fixed edge 
support conditions; T — 0 . 
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(a) Uniform temperature profile heating. 



(b) Dome-shaped temperature profile heating. 



(c) Roof-shaped temperature profile heating. 


Figure 11. Buckled shapes of square panel ( a/b = 1) under different temperature profile heating; 4C fixed edge 
support conditions; T = 0. 
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Figure 12. Plots of buckling temperatures (T 0 ) cr as functions of boundary heat sink temperature TJT a for different 
temperature profile heating; fixed cases. 
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Figure 13. Plots of buckling temperatures (T 0 ) cr as functions of boundary heat sink temperature TJT a for different 
temperature profile heating; free cases. 
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Figure 14. Plots of buckling temperature magnification factors of the first and the second kinds { 77 , §} as 
functions of boundary heat sink temperature TJT 0 . 
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(b) Roof-shaped temperature profile heating. 


Figure. 15. In-plane deformed shape of square panel under different temperature profile heating; four edges fixed; 
u(x,y ) ^ 0 , v(x,y) * 0 at panel interior; half panel plot. 
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Figure 16. Distribution of x-displacement u(x,y ) in square panel (a/b =1) under dome-shaped temperature 
profile heating; T = 1 °F; 4S fixed; u(x,y) * 0 , v(x,y) * 0 at panel interior. 
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Figure 17. Distribution of y-displaccmcnt v(x, y) in square panel ( a/b = 1) under dome-shaped temperature 
profile heating; T Q = 1 °F; 4S fixed; u(x,y) * 0, v(x,y) * 0 at panel interior. 
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Figure 18. Distribution of x-displacement u(x, y ) in square panel (a/b = 1) under roof-shaped temperature profile 
heating; T Q = 1 °F; 4S fixed; u(x,y ) ^ 0, v(x,y) * 0 at panel interior. 
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Figure 19. Distribution of y-displacemcnt v(x, y) in square panel ( a/b = 1) under roof-shaped temperature profile 
heating; T n - 1 °F; 4S fixed; u(x,y ) ^ 0, v(x,y) * 0 at panel interior. 
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Figure 20. Distribution of o x (x,y ) in square panel ( «/ b = 1) under dome-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u(x,y) * 0, v(x,y ) * 0 at panel interior. 
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Figure 21. Distribution of o y (x,y) in square panel (a/b = 1) under dome-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u{x,y) * 0, v(x, y) * 0 at panel interior. 
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Figure 22. Distribution of r Xy (x,y) in square panel (a/b = 1) under dome-shaped temperature profile heating; 
T = 1 °F; 4S fixed; u(x,y) * 0, v(x,y) * 0 at panel interior. 



Figure 23. Distribution of o x (x,y ) in square panel (a/ b - 1) under roof-shaped temperature profile heating; 
T o = 1 °F; 4S fixed; u(x,y) ^ 0, v(x,y) * 0 at panel interior. 
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Figure 24. Distribution of o y (x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u(x,y) * 0, v(x,y) * 0 at panel interior. 
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Figure 25. Distribution of r Xy (x,y) in square panel ( a/b = 1) under roof-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u{x,y) ^ 0, v(x,y) * 0 at panel interior. 
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Figure 26. Distribution of a x (x,y) in square panel ( a/b = 1) under dome-shaped temperature profile heating; 
T o = 1 °F; 4S fixed; u(x,y) = v(x,y) = 0 everywhere. 



Figure 27. Distribution of o y (x,y) in square panel ( a/b = 1) under dome-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u(x,y) = v(x,y) = 0 everywhere. 
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Figure 28. Distribution of o (x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
T Q = 1 °F; 4S fixed; u(x,y) = v(x,y) = 0 everywhere. 



Figure 29. Distribution of o y (x,y) in square panel (a/b = 1) under roof-shaped temperature profile heating; 
T = 1 °F; 4S fixed; u(x,y) = v(x,y) = 0 everywhere. 
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